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ABSTRACT 


This  thesis  examines  the  range  error  propagation  and 
uncertainties  associated  with  bistatic  sonar  operations.  An 
equation  for  determining  range  to  the  target  from  the  receiver 
is  explored  for  feasibility  of  practical  applications.  This 
particular  equation  does  not  require  the  separation  distance 
between  the  source  and  receiver  however,  it  does  require  an 
assvimption  of  the  mean  sound  speed  over  the  two  paths  of  the 
signal  even  though  it  could  change  drastically  over  a  few 
miles.  This  thesis  explores  the  contribution  of  eccentricity 
on  the  bistatic  ellipse  and  the  associated  error.  Exeuuples 
demonstrate  these  effects  by  comparing  cases  of  unequal  mean 
sound  speeds  over  the  different  paths  at  different  values  of 
eccentricity . 
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I.  INTRODUCTIOM 

A.  BACKGROUND 

Early  ASW  relied  solely  on  active  sonar,  or  Echo  Ranging. 
A  pulse  of  sound  was  projected  into  the  ocean  and  if  a 
submarine  was  in  the  path  of  the  signal,  an  echo  would  return. 
The  round-trip  time  of  travel  of  this  signal  was  used  with  the 
sound  speed  in  the  ocean  to  determine  the  range  to  the  target. 
This  method  of  ASW  was  very  effective  until  submarines  began 
to  exploit  the  water  conditions  of  the  ocean  and  became 
evasive,  i.e.,  hiding  in  shadow  zones,  fronts  and  eddies,  etc. 
The  projected  pulse  from  the  active  platform  usually  could  be 
detected  by  the  submarine  at  a  substantially  greater  range 
than  that  at  which  the  active  platform  could  detect  the 
reflected  pulse  from  the  target  because  of  the  greater 
transmission  loss  (TL)  over  the  round  trip  propagation  path 
[Ref.  1}.  This  was  one  of  the  factors  leading  to  interest  in 
passive  sonar.  Additionally,  passive  sonar  often  utilized 
lower  frequencies  having  smaller  losses,  which  could  allow 
longer  detection  ranges  [Ref.  2]. 

The  structure  of  the  speed-of-sound  profile  leads  to  a 
deep  sound  channel  which  can  be  exploited  in  convergence  zone 
(CZ)  propagation.  This  propagation  mode  with  its 
intrinsically  lower  acoustic  losses  provides  longer  standoff 
ranges,  keeping  the  detector  out  of  torpedo  ranges  of  the 
submarine.  On  the  other  hand,  a  disadvantage  of  passive  sonar 

'  -  1  ■ 


is  that  it  requires  a  relatively  time-conswjning  process  of 
bearing  rate  analysis  known  as  Target  Motion  Analysis  (TMA) 
to  determine  range  to  a  target  [Ref.  3].  Because  of  this, 
active  sonar  operations  are  still  important,  particularly  for 
last  minute,  close-in  prosecutions  where  exact  range  must  be 
determined  quickly  for  weapon  deployment. 

As  technology  advanced,  passive  sonar  became  more 
sensitive,  permitting  longer  detection  ranges.  Submarines, 
however,  reduced  their  radiated  noise  levels  nullifying  the 
advantages  of  the  more  sensitive  sensors.  This  progressive 
technology  also  made  way  for  the  continued  improvement  of 
active  sonar,  extending  capabilities  to  much  lower 
frequencies,  thus  permitting  active  sonar  exploitation  of  the 
deep  sound  channel.  Bistatic  Active  Sonar  (BAS)  has  brought 
about  longer  range  detection  by  active  sonar,  and  with  it  the 
advantage  of  immediate  determination  of  the  range  to  the 

i 

target.  Consequently,  the  vulnerability  of  active  sonar 
platforms  has  been  reduced  with  the  technique  of  histatic 
sonar  operations. 

B.  BISTATIC  SONAR 

Conventional  active  sonar  is  monostatic,  meaning  tliat  the 
source  and  receiver  are  usually  on  the  same  platform,  as 
displayed  in  Figure  1.  A  pulse  is  generated  by  the  source, 
transmitted  to  the  target,  reflected  from  the  target,  and  the 
resultant  echo  transmitted  back  to  the  receiver.  Dire(;tional 
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Vigur*  1  Monostatic  Arrangement 
sonar  sensors  provide  an  estimate  of  the  bearing  of  the  target 
and  the  range  to  the  target  is  determined  indirectly  by  the 
equation 

Range  «  (1»1) 

Where  V  is  the  sound  speed  in  the  ocean,  and  t  is  the  round 
trip  travel  time  of  the  pulse  and  echo. 

Blstatic  sonar  has  a  separated  source  and  received  on  two 
different  platforms.  Figure  2  is  a  representative  bistatic 
arrangement  [Ref.  4].  The  receiver  is  positioned  away  from 
the  source  to  receive  the  reflected  signal  from  the  target. 
In  the  bistatic  case,  the  source,  receiver,  and  target  form 
a  triangle  which  complicates  the  problem  compared  to  the 
monostatic  case  because  the  range  (R)  is  more  difficult  to 
determine.  The  source  is  presumably  detached  and  away  from 
the  threat,  reducing  its  own  vulnerability.  An  active 
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fiqwem  2  Bistatlc  Arrangement 
signal,  or  "ping",  originates  at  tlie  source,  travels  to  the 
target,  is  reflected  from  the  target,  and  then  travels  to  the 
receiver  as  shown  by  the  arrows.  The  t'^tal  travel  distance 
of  this  signal  is  X.  The  receiver  also  receives  a  direct 
signal  from  the  so\irce  over  the  path  S.  The  angle  (y)  is  the 
separation  angle  of  the  soxirce  and  the  target  me;.s\ircd  at  the 
receiver.  The  ultimate  goal  of  such  an  operation  is  to  obtain 
the  bearing  and  range  of  a  target  with  respect  to  the 
receiver.  Since  the  bearing  is  directly  obtainable  from  the 
directional  sonar  receiver,  determining  the  ranga  (R)  is  the 
problem  that  is  to  be  treated  here. 


4 


II.  STATEMEiro  OF  PROBLEM 

A.  RAZTGE  EQUATION 

An  equation  7.ust  be  derived  that  solves  for  range  (R)  to 
the  target  in  terms  of  known  quantities.  For  our  purposes,  we 
shall  assume  that  the  total  travel  distances  of  the  acoustic 
signal  over  the  two  paths  between  source  and  receiver,  the 
distances  X  and  S  in  Figure  2,  are  known.  Then,  given  the 
bearing  y,  we  have 

a  »  X-R  (2.1) 

and,  from  the  law  of  cosines, 

«  R*  +  5*  -  2RS'C0SY.  (2.2) 

Substitution  of  the  first  equation  into  the  second  allows  the 
elimination  of  a, 

(X-R)^  •  R}  *  -2RSC03y  (2.3) 

Expanding  the  term  on  the  left  gives 

X^ -2RX*R^  •  R^  *S^ -IRScoay  j  (2.4) 

Rearranging,  we  have 

2RX-2RScoay  -  (2.5) 

Solving  for  R  results  in  [Ref.  5] 

R  - 

2  (Jf-Scosy) 

The  above  equation  gives  R  in  terms  of  S,  X,  and  y. 


(2.6) 


The  time  of  travel  of  a  signal  from  the  hource  directly  to 
the  receiver  over  the  path  S  is  tj,  and  the  time  of  travel  for 
the  indirect  route  X  is  t^.  The  difference  between  t^  and  tj 


C  “  ^s‘ 


(2.7) 


These  time  definitions,  combined  with  the  mean  sound  speed 
over  the  propagation  path,  determine  the  distance  that  the 
signal  has  traveled.  Use  of 


X  =  Vt,  and  S  *  Vt^ 


(2.8) 


yields  several  different  equations  involving  different 
combinations  of  the  variables;  V,  S,  X,  t^,  t^,  t  and  y  [Ref. 
6].  One  such  combination  can  be  derived  as  follows. 
Substitution  of  (2.9)  into  (2.'6),  results  in 


(VC^)*-<V't.)» 

2  iVt^-Vt^r.osy)  ' 


(2.9) 


and,  the  substitution  of  (2.7)  for  t,,  provides 

2(V(t+C,)-V^C,COSY)  ’ 


(2.10) 


Expanding  both  the  numerator  and  the  denominator  leaves 


■*-2V“cc,  + 

2  (Vt-*- VC.- Vt.cosy  )  ' 


(2.11) 


and  simplification  gives 


6 


V*t*  +  2V*£t. 


2V(t*t,-  t^cosy )  ‘ 


Now,  use  of 


1  -  cosy  =  2sin^-|^ 


produces 


R  = 


Vt(Vt  +  2Vt,) 


2V(t+t,  -t,  +  2  t,sin^  ^  ^  ) 

A 


and  therefore 


R  = 


t(  Vt  +  2V't,) 


2  ( t  +  2  t,sin*y/2 ) 


Now,  use  of  the  relationship 


_£ 

tc 


gives 


IfiiSl 


2(t*2  t,sin*y/2 ) 


and  finally 


St(t*2t,) 


2  tg  (  c  ♦2  c^sin*y/2 ) 


(2.12) 


(2.13) 


(2.14) 


(2.15) 


(2.16) 


(2.17) 


(2.18) 


This  equation,  commonly  encountered,  is  a  function  of  S, 
t,  tj,  and  y.  Measured  quantities  of  the  variables  used  above 
are  often  known  to  much  more  precision  than  the  value  of  the 
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mean  sound  speed  V  over  many  miles.  Thus  V  has  the  largest 
amount  of  uncertainty,  because  the  sound  speed  profile  at  the 
receiver  must  be  assumed  over  the  whole  path  of  the  signal. 
(Most  bi static  sonar  operations  use  a  formula  similar  to 
Equation  2.18  to  determine  range  from  the  receiver  to  the 
target  [Ref.  6].)  However,  this  thesis  explores  the  use  of  an 
equation  that  contains  V  explicitly  rather  than  S  and  also 
allows  the  treatment  of  different  mean  sound  speeds.  By  using 
the  variable  S,  there  is  an  inherent  assumption  that  S  has  a 
readily  obtainable  value.  To  see  the  advantages  using  V 
rather  than  S,  it  is  appropriate  to  consider  the  role  of  ASW 
Command  Control  and  Communications  (C^) ,  few  operational  facts 
must  be  discussed. 

B.  CONSIDERATIONS 

Since  S  is  the  distance  between  the  source  and  receiver, 
an  accurate  position  of  the  source  must  be  known  for  the 
calculation  of  the  separation  distance,  S.  This  requires  that 
reliable  communications  be  maintained  between  the  two  vessels 
carrying  source  or  receiver  (Figure  3) .  anips  that  are  within 
"line  of  sight"  usually  communicate  by  tJHF  (Ultra  High 
Frequency)  radio,  which  is  extremely  reliable,  but  veJry 
limited  in  range.  When  ships  are  separated  by  a  distance 
greater  than  "line  of  sight",  approximately  20  nm,  other  means 
of  communication  must  be  utilized.  One  such  method  is  the  HF 
(High  Frequency)  radio  which  has  a  longer  range  because  of  the 
radio  wave  propagation  that  follows  the  curvature  of  the  earth 


Figure  3  Communications  Illustration 


[Ref.  7].  The  disadvantage  of  HF  communications  is  that  this 
type  of  communication  is  extremely  dependent  on  atmospheric 
conditions,  since  atmospheric  "ducting"  is  required.  Another 
disadvantage  of  HF  communications  is  the  vulnerability  of  the 
HF  transmitter  to  direction  finding  (DF)  [Ref.  8]  by  an 
opponent.  For  example,  there  are  Soviet  HFDF  sites  around  the 
world  that  can  triangulate  HF  transmissions  ro  determine  the 
location  of  the  transmitter.  A  third  method  of  communicating 
over  great  distances  is  UHF  Satellite  communications  (SATCOM) . 
In  this  method,  the  UHF  transmission  is  "bounced"  off  of  a 
satellite  to  another  station  over  the  horizon.  SATCOM  is 
extremely  reliable,  and  is  indeed  a  preferred  method,  but  in 


the  age  of  “'killer  satellites,"  dependency  on  satellites  can 
be  extreicely  limiting  in  operations. 

If  (2.18)  were  being  used,  then  a  communication  link  must 
be  maintained.  If  this  link  is  broken  and  communications  arc 
lost,  there  are  two  options  to  continue  bistatic  operations. 
The  first  option  is  that  the  source  and  receiver  move  to 
within  line-of-sight  UHF  range,  typically  20  nm.  This  would 
reduce  the  advantage  of  the  bistatic  concept  because  the 
receiver  is  not  "out  front"  exploiting  the  submarines 
reflected  energy,  and  decreasing  the  source's  vulnerability. 
The  second  option  in  case  of  lost  communications  is  for  the 
receiver  ship  to  approximate  the  source's  position  by  applying 
the  last  known  course  and  speed  to  the  last  known  position, 
also  known  as  Dead  Reckoning  (DR) .  This  introduces 
significant  error  in  the  variable  S. 

An  alternate  means  of  locating  the  target  using  bistatic 
operations  would  be  based  on  an  equation  that  removes  the 
necessity  to  )cnow  the  distance  S  directly. 

C.  AN  ALTERNATIVE  EQUATION 

To  avoid  the  difficulties  mentioned  in  the  previous 
section,  we  return  to  (2.6), 

- - 

2  (Jf-5cosy) 

and  use  the  relationships 

X  •  Utjf  and  S  •  vtg  (2.19) 
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to  introduce  differences  in  the  mean  sound  speed  over  the 
different  paths,  where  U  is  the  siean  sound  speed  over  the 
reflected  path  and  V  is  the  mean  sound  speed  over  the  direct 
path.  Substituting  (2.19)  into  (2.6)  results  in 


(C7t,)*- 

2  (Ut^-VtgCosy)  ' 


(2.20) 


and  by  dividing  the  numerator  and  denominator  by  it 

becomes 


4(1 


(2.21) 


A  scaled  time,  can  be  introduced  to  deal  with  these 

differences  in  speed,  where 


(2.22) 


Euqation  (2.21)  can  now  be  rewritten  as 


R  - 


tl) 

2(t,-t.cosY) 


(2.23) 


The  equation  still  requires  values  for  r^,  t,,  and  y,  but 
S  is  no  longer  directly  used.  This  thesis  examines  the 
utility  of  this  equation  (2.23)  as  an  alternate  to  (2.18). 

An  important  point  to  note  is  that  in  order  for  bistatic 
operations  to  take  place,  convergence  zone  conditions  must 
exist,  and  for  these  to  exist  there  must  be  a  minimum  of  200 
fathoms  of  depth  excess  for  reliable  propagation  [Ref.  3].  In 
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essence,  bistatic  operations  usually  must  take  place  in  the 
deep  ocean  and  away  from  typical  fronts  and  eddies  where  water 
conditions  change  drastically  over  short  distances. 

The  feasibility  of  using  this  equation  will  be  examined  by 
means  of  error  analysis.  We  include  in  this  analysis  the 
variation  in  mean  sound  speed  between  the  two  different  paths 
as  well  as  the  uncertainties  in  bearing  and  in  times  of  flight 
of  signal  and  echo. 


III.  DEVELOPMEMT  07  AMMiYSIS  MODEL 


Analysis  of  (2.23), 

j;- 

2  (Tjj- t^cosy)  ' 

is  conducted  from  the  point  of  view  that  it  is  the  equation  of 
an  ellipse.  Indeed,  any  bistatic  sonar  geometry  can  be 
described  in  terms  of  ellipses,  the  properties  of  which  relate 
directly  to  the  sonar  parameters. 


A.  ZMTRODnCTIOM  OF  ELLIPSE 

Figure  4  is  a  diagram  of  a  bistatic  arrangement  with  a 
superimposed  ellipse,  where  •  is  half  the  length  of  the  major 
axis  (the  semi^major  axis),  and  S  is  the  distance  between  the 
two  foci.  The  eccentricity  (e)  of  an  ellipse  can  be  defined  by 
the  ratio  e»S/2a.  Use  of  the  previously-defined  relationships 
and  properties  of  ellipses  then  provide  additional  es^ressions 


S  ,  5  ,  .  _£5 

X  a*R  Utj  T,  ■ 


(3.1) 


Equation  (3.1)  reveals  that  the  effect  of  a  difference  between 
the  mean  sound  speeds  U  and  V  over  the  two  paths  can  be 
completely  accounted  for  by  making  an  adjustment  in  the 
eccentricity  e  of  the  ellipse.  The  manipulations  ^at  follow 
will  result  in  convenient  expressions  for  the  fractional 
errors  in  range  arising  from  uncertainties  in  timing  errors , 
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differences  between  U  and  V,  and  bearing  errors.  The  value  of 
eccentricity  may  range  between  0  and  1,  where  e  «  0  is  a 
perfect  circle  and  e  -  1  is  a  very  flat  ellipse  approaching  a 
straight  line  joining  the  foci.  These  elliptical 
relationships  will  be  useful  in  describing  the  error 
associated  with  different  geometries. 


rigiira  4  Elliptical  Arrangement 
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B.  DIFFEREMTZATING  TERMS 

The  error  in  range,  AR,  can  be  expressed  in  terms  of  the 
errors  in  V,  tj,  and  y  as 


AR 


dv 


3J?  a  ^  ^  A  m* 


dt 


ay 


(3.2) 


or,  in  the  form  of  a  fractional  error,  as 


~R 


Rdv  .  R9t, 


At  ^ 


1  dR 

R  at. 


At.4fAr. 


(3.3) 


The  coefficients  of  the  error  increments  AV,  Ar^,  At,,  and  Ay 
on  the  right  side  of  the  equation  may  be  obtained  by 
logarithmic  partial  differentiation  of  (2.23).  The 
coefficient  of  AV  in  the  first  term  on  the  right  is 


IM  m  1 

R  dv  V' 


(3.4) 


The  coefficient  in  the  second  term  is 

IJI  M  _  1 

x^-tjcosy' 


(3.5) 


The  coefficient  in  the  third  term  is 

1  dR  ‘Zt.  ^  cosy 
Rdc,  T,»-t,»  T.-t.cosy' 


Finally  the  coefficient  in  the  y  term  is 


(3.6) 
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(3.7) 


C.  MEAN  SQUARE  ERROR 

The  fractional  error  in  range  can  be  obtained  through  the 
use  of  this  equation  if  all  errors  in  the  variables  V,  Tj,  tj, 
and  y  are  )cnown  completely  ,  both  in  magnitude  and  in 
direction.  On  the  other  hand,  if  all  errors  (i.e., 
uncertainties)  in  the  variables  V,  t,,  and  y  are  assumed 
to  be  statistically  independent  of  each  other,  they  may  be 
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combined  into  a  mean  square  error  (KSE) .  The  HSE  is  used  to 
determine  the  total  combined  error  resulting  from  errors  in 
variables  in  the  range  equation.  It  is  calculated  by  taking 
the  square  root  of  the  sum  of  the  squares  of  the  individual 
errors.  Each  of  the  terms  has  an  Uncertainty  value  (A) 
associated  with  that  particular  variedsle,  i.e.,  At^.  The 
correct  numerical  uncertainty  values  of  each  variable  used 
operationally  has  a  security  classification  and  will  not  be 
used  in  this  thesis.  To  simplify  the  notation  in  the  model 
the  coefficients  of  the  error  terms  in  (3.8)  will  be  expressed 
as  Ay,  Aj  or  Aq,  where 


Ax 

Ag 

Ag 


(.  2  .  _  1  ] 
Vl-e’  l-ecosyj' 

V l-ecosy/ 

~e  slny 
1-ecosy  ’ 


The  use  of  these  coefficients  gives  the  fractional  error  in 
range  in  the  form 


AR 

R 


(3.10) 


The  mean  square  fractional  error  is  thus 
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1(A«, 


2.  (3.11) 


ir  which  it  is  assumed  that  the  errors  on  the  right  side  of 
the  equation  are  statistically  independent.  Equation  3.11 
calculates  error  in  terms  of  the  fractional  timing  errors 
referenced  to  a  time  base  of  t,(.  If  it  is  preferable  to 
express  timing  errors  over, each  path  referenced  to  the  travel 
time  over  each  path,  then 


AR  ^ 


+  eAt 


Ac« 


(3.12) 


Not0l 


A  s 

Vtg 

1 

=  i?eJL 

Ut/ 

V  ts 

^A 

.  ^A 

t. 

-  eAg 


A  Cg 


The  mean  square  fractional  error  is  then  changed  to 


N 


♦V(Ay)» 


(3.13) 


To  understand  the  use  of  these  equations  suppose  we 
consider  the  general  method  by  which  the  range  from  the 
receiver  ship  to  the  target  (R)  and  the  associated  range  error 
can  be  found.  Assume,  for  example,  that  the  signal  travel 
times  ty^  and  t,,  and  the  azimuth  angle  y  all  have  been 
measured.  Assume,  further,  that  the  average  sound  speed  V 
over  the  direct  path  from  source  to  receiver  is  known.  Also 


18 


let  U,  the  average  sound  speed  over  the  secondary  path, 
source-to-target-to-receiver,  be  the  same  as  V,  so  that 
The  range  R  can  then  be  calculated  from  2.23,  where 

^  ^  V(tl-tl) 

2  (tx-  t^cosy) 


Assume  also  that  the  uncertainties  in  time,  At^  and  At^, 
and  the  uncertainties  in  azimuths.  Ay,  are  all  known.  The 
fractional  error  (uncertainty)  in  range  associated  with  each 
of  the  uncertainties  in  the  separate  variables  can  be  obta'<ned 
as  follows. 

For  timing  error  in  t,  alone  the  corresponding  fractional 
error  in  range  is 


(3.14) 


He  obtain  A,  from  the  graphs  of  Fig.  5  through  13  after  first 
calculating  the  eccentricity. 


for  this  particular  case  being  studied.  Similarly  for  error 
in  t,  alone  use 


/  AR 

\  R 


). 


At, 


OAm 


At, 


(3.15) 


and  for  error  in  y  use 
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where  Ay  must  be  in  units  of  radians.  If  more  than  one  of  the 
independent  variables  is  uncertain,  then  either  (3.11)  or 
(3.13)  for  the  combined  errors  should  be  used. 

D.  ERROR  ANALYSIS 

A  computer  program  was  written  to  study  the  behavior  in 
range  error  in  each  of  the  three  variables  t,,  tj,  and  y.  The 
program,  written  in  FORTRAN,  was  used  for  nine  different 
values  of  eccentricity  a  of  the  basic  range  ellipse,  from  0.1 
to  0.9.  The  output  data  of  the  program,  the  values  of  A^,  Aj, 
and  Aq,  were  plotted  in  the  fora  of  curves  from  which  error 
estimates  can  be  made  for  any  reasonable  scenarios. 


ZV.  RESULTS  AMD  VERIFICATION 

A.  RESULTS  OF  MODEL 

Figures  5  through  13  are  plots  of  the  output  data  of  the 
previously  mentioned  Fortran  progreun.  All  of  these  plots 
Indicate  some  relationships  which  will  be  examined  further. 
Each  of  the  nine  figures  is  a  graph  of  the  error  coefficients 
for  a  different  eccentricity  m,  beginning  with  Figure  5  for 
the  case  of  eccentricity  0.1  and  continuing  to  Figure  13  with 
eccentricity  of  0.9.  In  all  of  the  plots  the  values  are 
represented  by  solid  lines,  the  A^  values  are  represented  by 
dotted  lines  and  the  Aq  values  are  indicated  by  dashed  lines. 
The  error  coefficients  (ordinates)  give  the  magnitudes  of  the 
errors  associated  with  each  of  the  three  terms.  The  baseline 
number  (S)  used  for  input  is  inconsequential  because  the  value 
of  a  chafacterizes  the  driver  in  the  error  propagation.  Each 
of  the  ttiree  terms  will  be  evaluated  Individually. 

The  tenti  is  the  indicator  of  the  behavior  of  fractional 
error  associated  only  with  uncertainty  in  the  value  of  t^. 
The  graplis  in  the  figures  show  how  A^  and  the  other  error 
coefficicsnts  vary  with  bearing  angle  y  The  graphs  reveal 
that  for  low  values  of  eccentricity  a,  the  error  coefficient 
begins  small  and  shows  a  small  increase  as  y  increases  to 
180*.  For  large  values  of  a,  the  values  of  A^  remain  small 
for  ibw  values  of  y,  but  as  y  increases,  the  A^  value 


21 


increases  rapidly  until  about  mid  range  of  y,  approximately 
70*,  when  it  levels  off  and  ap»5roaches  a  maximum  at  180*. 
This  information  indicates  that  for  small  values  of  e,  the 
error  contribution  by  the  term  is  small  for  all  values  of 
y.  On  the  other  hand,  for  large  values  of  eccentricity,  the 
error  contribution  of  t^  begins  small  for  angles  less  than 
10*,  but  rapidly  increases  for  larger  angles  until  it  steadies 
out.  In  practical  terms,  recall  that  a  small  value  of  e  means 
a  small  separation  distance  (S  ) between  the  source  and 
receiver  compared  to  the  distance  between  the  receiver  and  the 
target  (R) .  Consequently,  a  smaller  S  means  a  smaller  error 
associated  with  t^,.  For  larger  values  of  e,  the  error 
contributed  by  t^  increases  drastically  for  separation  angles 
(Y)  greater  than  10*. 

The  Kg  coefficient  characterizes  the  effect  of 
uncertainty  in  the  time  tg.  These  values  are  represented  by 
the  dotted  line  on  each  of  the  nine  graphs.  For  y  of  0*  and 
180*  this  error  coefficient  is  about  the  same  as  for  the 
previously  determined  Ax.  For  small  values  of  e,  Aj  gradually 
decreases  to  zero,  at  about  70*,  and  then  increases  until  it 
has  the  same  value  as  Ax  at  180*.  For  higher  values  of  e,  the 
value  of  Aj  drops  more  rapidly  to  zero  at  angles  less  than  10* 
and  then  approaches  the  values  of  Ax  as  it  increases  sharply 
and  behaves  asymptotically  done  about  70*.  The  Aj  coefficient 
goes  to  zero  at  moderate  values  of  y,  about  70*,  for  low  •, 
and  it  goes  to  zero  at  small  angles  for  large  e.  When  the 


22 


source  and  receiver  are  closer  to  each  other  than  the  receiver 
is  to  the  target,  this  angle  at  which  the  error  vanishes  is 
around  70*  to  80*,  and  for  large  values  of  e,  this  angle  is 
much  lower,  near  10*. 

The  error  coefficient  A^,  associated  with  the  y  term,  is 
zero  for  y  values  of  0*  and  180*,  as  expected.  For  smaller 
values  of  e,  the  error  coefficient  increases  slowly  until  it 
peaks  at  about  90*.  For  small  values  of  e,  the  contribution 
of  the  y  remains  small  compared  to  the  other  two  terms.  For 
larger  values  of  e,  the  maximum  .  value  of  A.  increases 
significantly,  arid  it  "shifts'*  to  smaller  angles  to  the  left 
as  e  increases.  In  practical  terms,  when  the  source  is  closer 
to  the  receiver  than  the  target  is  to  the  receiver,  that  is, 
when  S  is  less  than  R,  the  error  contribution  of  the  y  value 
is  relatively  insignificant.  If  the  geometry  is  such  that  • 
is  large,  then  the  error  would  be  reduced  if  the  angle  y 
remains  greater  than  90*.  All  plots  shown  in  the  figures 
represent  the  absolute  values  of  the  error  coefficients. 
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rigura  6  Error  vs.  Angle  for  e«0.2 
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B.  VERIFICATION  OF  MODEL 

Verification  of  this  model  was  done  by  graphical 
confirmation  (elliptical)  of  theoretical  expressions.  The 
three  varizdiles  that  contribute  errors  to  a  particular 
bistatic  geometry  are  t,,  t^,  amd  y.  The  variadsles  t,  and  t, 
determines  the  eccentricity  of  the  ellipse  and  y  determines 
the  location  of  the  target  on  that  ellipse.  For  each  of  the 
three  curves ,  one  of  the  three  variables  is  assumed  to  be  in 
error  by  approximately  10%  while  the  remaining  two  variables 
were  considered  to  be  exactly  determined.  This  forces  the  two 
remaining  terms  to  become  zero,  and  the  fraction«.l  error, 
AR/R,  is  now  equal  to  the  one  remaining  term.  The  fractional 
error  can  be  crudely  determined  from  the  elliptical  figure  at 
various  values  of  y  and  compared  with  the  value  of  the 
appropriate  coefficient  A,,  A,,  or  Ag  obtained  from  one  of  the 
figures  (Figures  5-13)  for  a  corresponding  value  of  e.  The 
first  of  the  terms  to  be  validated  is  the  A,  term.  Two 
ellipses  are  drawn  in  Figure  14  that  share  the  same  focal 
points,  because  t,  is  held  constant,  thus  guaranteeing  a 
constant  inter-focal  distance  S.  Ellipse  1  has  an 
eccentricity  of  0.80.  If  the  t^  value  is  increased  by  10%, 
the  eccentricity  can  be  seen  from 

C- 

e  -  for  U"V  (4.1) 
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a  small  value  of  y,  the  AR/R  should  be  relatively  small  since 
R  has  its  largest  values  here.  As  y  Increases,  AR  increases 
as  R  decreases  resulting  in  a  rapid  climb  of  AR/R  which  is 
confirmed  in  Figures  11  and  12.  Finally,  as  y  approaches 
180*,  AR  and  R  remain  constant  which  explains  the  asymptotic 
behavior  of  the  function  at  higher  values  of  y.  , 

The  next  term,  Aj,  is  validated  through  Figure  15,  where 
t,  and  y  are  held  constant  and  tj  is  decreased  by  10%.  Ellipse 
1  in  Figure  15  is  the  original  geometry  for  a  given  tj  and 


Figure  15  Validation  of  A, 

ellipse  2  shows  the  geometry  with  t,  reduced  by  10%,  which 
gives  a  reduced  eccentricity.  This  graphic  is  moat 
interesting  because  there  is  an  immediate  explanation  for  the 


zero  value  of  A3  at  about  20*  in  Figure  11.  Shortly  beyond 
this  value  of  y,  the  plot  behaves  quite  similarly  to  the  A^ 
plot.  As  Y  increases  beyond  about  90*,  the  graph  again 
demonstrates  asymptotic  behavior. 

The  validation  of  the  A^  function  is  shown  in  figure  16. 


Figure  16  Validation  of  A^ 

The  validation  in  this  case  requires  that  t^  and  t,  remain 
constant  while  y  is  permitted  to  change  slightly.  Since  t^ 
and  tj  remain  constant,  so  must  e  remain  the  same.  Thus  there 
is  only  one  ellipse  to  be  considered  here.  Figure  16  shows 
four  angles  y  and  four  associated  Ly*a  shown  by  the  dotted 
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lines.  The  difference  in  len^h  between  the  dotted  line  and 
the  solid  line  is  the  AR  value.  Note  that  AR  is  zero  where  y 
is  either  0*  or  180*.  The  graphic  demonstrates  that  for 
higher  values  of  e,  as  in  this  case,  there  will  be  a  maximum 
value  of  AR  between  20*  and  50*,  followed  by  a  gradual  decline 
to  zero  as  y  increases  to  180*.  For  lower  values  of 
eccentricity,  the  AR  value  reaches  its  maximum  close  to  90* 
and  has  a  much  smaller  value  than  in  the  case  of  larger 
eccentricity. 

The  geometrical  analysis  resulting  in  Figures  14  through 
16  gives  a  graphical  interpretation  of  the  errors  arrived  at 
algebraically  in  Chapter  III.  From  these  figures  one  can 
readily  see  which  errors  are  more  significant  for  the  various 
conditions  that  may  occur  in  the  range  problem. 
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V.  EXAMPLE  CASES 


It  is  helpful  to  illustrate  use  of  the  analysis  by 
applying  it  to  practical  problems.  The  following  examples 
will  demonstrate  the  previously  mentioned  error  propagation  in 
cases  that  might  be  encountered  at  sea.  These  examples  will 
be  approached  as  if  the  calculation  were  conducted  from  input 
provided  by  Sonar  Control  onboard  an  ASW  ship.  The 
uncertainty  numbers  used  in  these  examples  are  not  necessarily 
realistic,  but  are  useful  in  the  demonstration  of  the  method. 
The  first  three  examples  are  of  geometries  of  low  eccentricity 
while  the  remaining  examples  are  of  high  eccentricity. 

A.  EXAMPLE  l:  LOW  ECCENTRICITY;  NO  UNCERTAINTIES  IN  V 

For  the  first  example  assume  that  the  mean  sound  speed  is 
the  same  over  both  paths,  that  is,  U  ecjuals  V.  Figure  17  is 
a  diagram  of  an  elliptical  arrangement  of  low  eccentricity;  S 
is  less  than  X.  Suppose  that  the  following  data  are  provided: 

Y  =  110*  ±1* 

42.62  ±  0.01  5 
Cg  «  213.9  t  0.03  s 
V  -  4990  fC/s 
C7  -  V 

conversion  ■  6076  ft/nm 

Calculated  values  of  the  various  terms  in  the  equations  are 
then  as  follows: 
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- ]  =  1-147 

Vl-e^  1-ecosYj 

A,  °°°1f  ]■  -0.737 

^l-e“  l-ecoey J 


-e  sxny  , 

“  1-ecosY 

If  it  is  asstuned  that  the  uncertainties  are  random 
(independent)  then  the  RMS  fractional  error  may  be  calculated 
from  (3.11)  in  which  it  is  assumed  that  U  =  V,  so  that  =  tj,, 
and  also  that  AV/V  equals  0: 

-|(A/?)^„-^  -"V (Ay) *  =  0.003075 

Thus  the  expected  uncertainty  in  range,  AR,  should  be  about  ± 
500  yards  or  less. 

B.  EXAMPLE  2:  LOW  ECCENTRICITY;  UNCERTAINTY  IN  V 

Assvime  the  same  situation  as  in  Example  1,  but  now  assume 
an  uncertainty  in  V  of  ±  3  ft/s.  Now,  we  have  an  AV/V  value 
of  0.0006  in  addition  to  the  other  errors  in  Exa^uple  1.  The 
term  AV/V  in  (3.11)  is  now  not  zero,  but  «  t^  as  before. 
The  fractional  error  is  now 

^  ( A J!) »  -  ^  (■^F  .V  (At)  '■  0 . 003133 

Note  that  the  error  in  AV/V  in  this  case  does  not  have  a  large 
effect. 
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C.  EXAMPLE  3t  LOW  ECCENTRICITY;  U  AND  V  UNEQUAL 

In  this  example  the  situation  is  the  seme  as  in  Example  1 
but  now  let  U  -  5040  ft/s  so  that  U  is  not  equal  to  V.  The 
data  to  be  used  are  in  this  case: 

Y  =  110*  ±1* 

=  42.62  ±  0.01  s 
Cjf  -  213.9  ±  0.03  s 

V  *  4990  ft/ 3 
U  »  5040  ft/ 3 

Nov,  it  is  necessary  to  rescale  t^  as  discussed  in  Chapter 
II,C. , (2.22) : 

T,  -  -1.010 

From  this  point,  analysis  proceeds  as  in  Example  1. 

e  -  ^  -  0.198 


2  {Xj-  tgcoay) 


79.873nm 


Ajf  -  (— - - - ]  -  1.145 

Vl“«*  1-ecosY^ 


-0.733 


(  -2 e  ^  COSY  \  . 
^l-e^  1-ecosy  j 


^0 


-e  siny 
1-ecosy 


-0.174 


Now,  we  use  (3.11)  with  r,,  *■  t,,  and  AV/V  «  0,  to  get 


^  ( AJ?)  ^  ( Ay) »  -  0 . 003133 

The  mean  square  fractional  error  is  now  somewhat  greater  than 
that  in  Example  1,  as  might  be  expected. 

D.  EXAMPLE  4:  HIGH  ECCEMTRICXTT;  MO  UMCERTAZMTZE8  ZM  V 

Assume  that  the  mean  sound  speed  is  the  same  over  both 
paths.  Figure  18  is  a  diagram  of  an  elliptical  arrangement  of 
high  eccentricity,  which  means  that  the  source  and  receiver 
are  generally  farther  apart  than  the  target  and  receiver. 
Suppose  that  the  following  data  are  provided: 


y  »  liO*  ±  1* 

>  85.23  ±  0.01  s 
tg  -  106.54  ±  0.03  s 
V  -  4990  ft/ 3 
U  •  V 
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^  o  *  12.37 mn 

2(t,-t5CosY) 


Target 


Pigure  18  Example  Case,  High  Eccentricity 

Aji  -  f— ^"1 - — — ]  “  4.770 

Vl”®*  1-ecosYj 


4.713 

VI"®*  l-®cosYj 


:5-8inL  .  -0.590 


"  1-ecosY 
Equation  (3.11)  gives  for  this  case: 


(Ay) *-0.010398 

Note  that  AR/R  is  almost  a  factor  of  3  greater  than  in  the 
previous  examples  of  low  eccentricity.  However,  for  the 
geometry,  the  actual  range  error  is  260  yards,  about  half  of 
the  expected  error  of  the  earlier  cases.  It  is  Important  to 
point  out  that  the  fractional  error  AR/R  is  a  fairly  sensitive 
fuctlon  of  eccentricity,  but  the  actual  range  error  depends  on 
the  size  of  the  ellipse. 
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1.  EXAMPLE  5:  UMCERTAIMTY  IM  ▼ 

Here  the  same  data  are  provided  as  in  Example  4,  but  in 
this  case  V  has  uncertainty  of  ±  3  ft/s.  This  time  we  have 
AV/V  s  0.0006  to  be  included  in  the  RMS  calulation,  as  was 
used  in  Example  2,  and  again  =  t^.  The  RMS  fractional 
uncertainty  is  thus 


The  expected  range  error  in  this  case  is  also  about  260  yards. 


F.  EXAMPLE  «:  U  AMD  V  UNEQUAL 

In  this  example,  as  in  Example  3,  assume  that  the  same 
data  are  provided  as  in  the  two  previous  examples  except  that 
U  is  not  equal  to  V  and  there  is  no  uncertainty  in  V.  Thus, 

Y  *  110*  ±  1* 

tg  =  85.23  ±  0.01  s 
tg  -  106.54  ±  0.03  3 

V  •  4990  ft/ 3 
U  •  5040  ft/ 3 

Again  we  must  rescale  t^  by  using  r^: 

Xg  ■  -^tg  ■  1.  Olt, 

and  the  process  proceeds  as  in  Example  4. 

e  -  -  0.792 

Use  (3.11)  with  tg  *  Tg,  and  AV/V  ■  0  to  get  a  fractional  error 
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-0.010327 

The  expected  range  error  is  here  about  270  yards,  slightly 
larger  than  in  Examples  4  and  5. 
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VI .  CONCLUSION 


A.  DISCUSSION  AND  RESULTS 

Improved  quieting  techniques  by  submarines  have  been  very 
effective  in  reducing  the  exploitable  emissions  for  passive 
sensors.  On  the  other  hand,  detection  by  active  sonar  is 
independent  of  the  submarine's  quieting  efforts  and  gives  the 
added  bonus  of  providing  an  immediate  determination  of  range 
to  the  target.  Thus,  active  sonar  seems  to  give  certain 
obvious  advantages.  Additionally,  a  bistatic  sonar  technique 
increases  detection  range  of  active  sonar  and  reduces  the 
vulnerability  of  the  active  platform. 

This  thesis  has  examined  a  method  of  predicting  range  and 
range  uncertainties  in  bistatic  sonar  operations.  The  method 
described  here  offers  flexibility  over  conventional  methods 
that  rely  heavily  on  communication  between  sovurce  and  receiver 
for  successful  operations.  Operational  methods  that  use  the 
separation  distance  S  are  often  preferred  over  those  that 
require  a  knowledge  of  the  mean  sound  speed  V,  especially  near 
fronts  and  eddies  where  the  mean  sound  speed  may  change 
drastically  in  a  few  miles.  However,  this  thesis  has  shown 
that  the  preference  of  using  the  baseline  over  the  mean  sound 
speed  is  unfounded. 

As  shown  in  the  previous  chapter,  the  range  error  is 
controlled  to  a  great  extent  by  the  eccentricity  of  the  range 
ellipse,  the  largest  errors  occurring  in  the  cases  of  high 


eccentrlcicy .  Examples  are  given  to  provide  representative 
applications  of  the  range  Equation  (2.3)  and  to  give  the 
estimates  of  the  range  errors. 

We  have  seen  that  fractional  error  AR/R  tends  to  be  small 
when  •  is  small,  although  dependance  on  •;  is  still  quite 
significant.  This  suggests  that  small  base-line  geometries 
(small  values  of  S  leading  to  small  •)  are  to  be  preferred  if 
these  are  operationally  feasible. 
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